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I. INTRODUCTION 

The derivation of low-energy hadronic observables, e.g. 
meson masses and decay constants, from the theory of 
the strong interaction (QCD) has so far proven to be im- 
possible by means of analytical methods. Because of this 
situation, numerical Lattice QCD simulations, whereby 
the functional integral is evaluated numerically on a dis- 
cretized space-time lattice, have developed into a major 
field of study. Such simulations are, however, seriously 
hampered by difficulties in the simulation of dynamical 
sea quark effects. Although much progress has been made 
recently, it is still impractical, for computational reasons, 
to simulate with sea quark masses that are close to the 
physical u, d quark masses of a few MeV. This situation, 
with sea quark masses of a few tens of MeV, is referred 
to as partially quenched (PQ) QCD. Consequently, the 
physical values of the sea quark masses have to be reached 
by extrapolation from the partially quenched simulation 
results. 

A practical method for this extrapolation is provided 
by Chiral Perturbation Theory (xPT), which provides 
the correct quark mass dependences of the various physi- 
cal quantities that are measured on the lattice. Standard 
three-flavor xPT as introduced by Weinberg, Gasser and 
Leutwyler in Refs. JJ, is valid in the (unquenched) QCD 
case of equal valence and sea quark masses. The general- 
ization of xPT to the quenched case (without sea quarks) 
or to the partially quenched case (sea quark masses dif- 
ferent from the valence ones) has been carried out by 
Bernard and Golterman in Refs. 1113. The quark mass 
dependence of partially quenched chiral perturbation the- 
ory (PQxPT) is explicit, and thus the limit where the sea 
quark masses become equal to the valence quark masses 
can be taken. As a consequence, is included in 

PQxPT and the free parameters, or low-energy constants 
(LEC:s), of yPT can be directly obtained from those of 

PQxPT II, 3. 

The calculation of charged pseudoscalar meson masses 
and decay constants to one loop (NLO) in PQxPT has 
been carried out in Refs. 0,0,0; and first results for the 
mass of a charged pseudoscalar meson at two loops or 
next-to-next-to-leading order (NNLO) in PQxPT have 



already appeared, for degenerate sea quark masses, in 
Ref. [3 . The need for such calculations is clear as NNLO 
effects have already been detected in Lattice QCD sim- 
ulations 01 111 • A calculation of the pseudoscalar meson 
masses for nondegenerate sea quarks is in progress 0. 

This paper presents the first calculation of the decay 
constants of the charged, or flavor off-diagonal, pseu- 
doscalar mesons in NNLO PQxPT, for three flavors of 
sea quarks (risea = 3). The results are characterized by 
the number of nondegenerate valence and sea quarks, de- 
noted dvai and dsoa, respectively. For the decay constants 
of the charged pseudoscalar mesons, the maximum num- 
ber of nondegenerate valence quark masses is dvai — 2. 
The degree of quark mass degeneracy in each result is 
sometimes also referred to with the notation c?vai + c^soa- 
The decay constant of the charged pion in the SU (2) sym- 
metric limit thus corresponds to the 1-1-2 case. Likewise, 
the decay constants of the charged and neutral kaons 
may be obtained from the dvai = 2 results with dgoa — 
2. Results are also presented for the case of dsea = 1 
(all sea quark masses equal), and dsea = 3 (all sea quark 
masses different). An extension of the present work to 
the neutral pseudoscalar mesons is also planned. 

The analytical expressions for the NNLO shift of the 
decay constants are in general very long, but the expres- 
sions simplify considerably when pairs of sea or valence 
quark masses become degenerate. In view of this, the 
NNLO loop results are given separately for each case of 
rfvai + c^sea Considered. In the next sections, the technical 
background for the NNLO calculations, the full results for 
the decay constants of the charged pseudoscalar mesons 
and numerical results as a function of the input quark 
masses are given, along with a concluding discussion. 



II. TECHNICAL OVERVIEW 

Most of the technical aspects that concern the calcu- 
lation of the pseudoscalar meson decay constants to two 
loops, or NNLO, are identical to those of the calcula- 
tion of the pseudoscalar meson mass, and have already 
been justified in Ref. 0. Most significantly, the La- 
grangians of PQxPT at 0{p^) and 0{p^) may be directly 
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obtained from the corresponding Lagrangians of normal 
(unquenched) np flavor xPT, provided that the traces 
and meson matrices are replaced with the supertraces 
and meson matrices relevant to the partially quenched 
theory 0, H H- This can be argued from the Replica 
method as in Ref. oi' by the fact that all the rela- 
tions used to constrain the number of terms in Ref. [TH 
remain valid when traces are changed to supertraces. We 
work here in the version of PQxPT without the <f>o as 
discussed in Ref. (5(1 . 



A. Masses and Low-Energy Constants 

All calculations in this paper have been performed with 
three flavors of valence quarks, three flavors of sea quarks 
and three flavors of bosonic 'ghost' quarks. These may 
be viewed as the u, d and s quarks in the valence, sea 
and ghost sectors, respectively. The purpose of the ghost 
quarks is to remove the valence quark loops which are 
disconnected from the external legs. 

The input quark masses rrii enter into the calculation 
in terms of the lowest order squared meson masses Xi, 
which are defined as usual in xPT, by Xi = '^Borrii. 
In the present calculations, we thus have three valence 
inputs XiiX2,X3, three sea inputs X4,X5,X6, and three 
ghost inputs X7: Xs, X9- In order for the disconnected va- 
lence quark loops to be canceled, the masses of the ghost 
quarks are always equal to those of the corresponding va- 
lence quarks, such that X7 = Xi, Xs = X2 and X9 = Xa- 
ExpHcitly, for dvai = 1, we have = ^2 = X3, for 
rival = 2, we have Xi = X2 7^ X3 and for dvai = 3, we have 
Xi 7^ X2 7^ X3- Similarly, for the sea quarks dsca. = 1 im- 
plies X4 = X5 = X6, while 4ca = 2 implies X4 = X5 7^ Xe 
and finally 4ea = 3 that X4 7^ X5 ^ Xe- 

The number of independent low-energy constants in 
unquenched and partially quenched xPT is slightly differ- 
ent, but the former are always linear combinations of the 
latter. For PQxPT, they are Fq and Bq at leading order, 
LI through L^^ at NLO and through Kl^^ at NNLO. 
In contrast, for three flavor unquenched xPT they are 

r(3) r(3) 

Fq and Bq at leading order, Lq^ through at NLO 
and C[ through Cg^ at NNLO. Note that the parameters 
Lli and correspond to the usual Fl^^- Also, the 

parameters Fq and Bq are identical for unquenched and 
partially quenched xPT. At NLO, the relations between 
the low-energy constants are [l|, [U Ull 

ri-(3) _ rr , rr 
^2 — ^0 + ''^2' 

■* = — 2Lq -I- L3, 

r'-(3) _ ('[\ 
^4.. .12 — ^4. ..12 ; 

and the corresponding linear relations relevant for the 
NNLO parameters can be found in Ref. 



B. Meson Propagators in PQ^PT 

The calculation of the pseudoscalar meson mass in 
Ref. was only performed for dvai = ^soa = 1- However, 
the present calculation of the decay constants to NNLO 
is also concerned with the more general cases of dvai and 
dsea = 2,3. This leads to much more involved expressions 
because of the appearance of the residues of the neutral 
meson propagators. We recall here the results of Ref. 
which we have translated to Minkowski space from the 
Euclidean formalism used there. 

In Minkowski space, the propagator Gf^ of a charged, 
or flavor-off-diagonal meson with flavor structure qiqj in 
supersymmetric PQ^PT, is given by fsj 

The factor Xij is defined in terms of the squared masses 
as Xij = iXi +Xi)/2- The sign vector ej is defined as -1-1 
for the fermionic valence and sea quarks (j = 1, . . . , 6) 
and —1 for the bosonic ghost quarks (j = 7,8,9). The 
propagator of a flavor-neutral meson in supersymmetric 
PQxPT is more complicated, as it connects mesons of 
different flavor indices as well. The propagator for a me- 
son with quantum numbers QiQi to one with quantum 
numbers qjqj may be written in Minkowski space as 

G-{k) = G^(fc)^,^_^G«,(fc). (3) 

The nontrivial part G^^- of the neutral meson propaga- 
tor may be expressed, by means of partial fractioning, in 
terms of a sum of single and double poles in fc^ . For the 
most general case of dsca = 3, in terms of the single-pole 
residue R, the double-pole residue R'^ and the auxiliary 
residue R'^, the propagator G^ is 
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+ 
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(4) 
(5) 



— Xtt + k"^ — Xi-) + *e 



For uniformity of notation, the squared (lowest order) 
masses of the neutral pion and the eta meson in the sea 
quark sector have been denoted by Xti- and Xrp respec- 
tively. They are functions of the sea quark masses, and 
in the case of dsca = 3 they are given by the lowest order 
xPT result with — r\ mixing active j^. They can be 
obtained from the relations 

Xtt + X>) = 2 (X4 + X5 + Xe) , 

XttXt) = \ (X4X5 + X5X6 + X4X6) ■ (6) 
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For dsea = 2, the 7r° pole in G'^ disappears, and in that 
case the (lowest order) eta meson mass is trivially given in 
terms of the remaining sea quark masses as = 1/3 X4+ 
2/3 xg. In this case the index tt has been suppressed in 
the residue notation. Further, for dsca = I7 both the tt'^ 
and 7] poles disappear and consequently both the tt and 
77 indices have been suppressed. In this way we obtain a 
unique notation for the residues, from which the number 
of nondegenerate sea quarks is immediately apparent. 

The various residues R of the propagators of the neu- 
tral, or flavor-diagonal mesons appear in the results, 
and are one reason why the PQxPT expressions are much 
more involved than the xPT results of Ref . • The use 
of the propagators 10} and |(SJ) in the present form has 
the advantage of producing results in terms of standard 
loop integrals which can be treated with known methods. 
On the other hand, the various residues R of the flavor- 
neutral meson propagator fulfill a very large number of 
relations and the direct output from the calculations of 
the diagrams consequently produces a large number of 
redundant terms. This problem does not yet manifest 
itself at the one-loop level, but becomes troublesome at 
the two-loop level when the mass degeneracies in the sea 
and valence quark sectors are lifted. With a major effort, 
the end result can be simplified, in some cases it has been 
compressed by more than an order of magnitude. 



C. Notation for Propagator Residues 

The form of the single-pole residues i?*^,; and the 
double-pole residue Rf, which appear in Eqs. ^ and ||SJ|, 
depends on the degree of degeneracy in the sea quark 
masses, which in turn is indicated by the number of in- 
dices in the single-pole residue R. It is useful to define 
the more general quantities R^ ^ such that 



^ab — 

pz _ 
^abc 



^abcd 



Xa - 
Xa 



Xb, 
■ Xb 



R: 



abode f g 



Xa - Xc 

{Xa - Xb){Xa - Xc) 

Xa~Xd 
iXa - Xb){Xa - Xc){Xa " Xd) 
iXa - Xe){Xa " Xf)iXa " Xs)' 



(7) 



and so on. The R^ notation is primarily useful for defin- 
ing the residues of the flavor-neutral propagators, but it 
may also appear independently in the final result for the 
decay constant. In such cases, the R^ have been gen- 
erated by simplification procedures, as all the residues 
that are naturally generated by partial fractioning of the 
propagator Gf^ are of the form given below. Note that 
R^ ^ has the same dimension as Xi for even num- 
ber of indices and is dimensionless for an odd number of 
indices. 

For the case of dsea = Ij all residues associated with 
the sea quark sector have reduced to numbers. Some 



nontrivial residues still appear if the valence quarks are 
nondegenerate, according to 



Rid ■ 



(8) 



For convenience, Rf is also used for dvai = 1 in order to 
maintain a notation similar to the more general cases. 

As noted in Ref. [3, the residues also simplify for 
dsca — 2. In particular, the residue of the neutral pion 
pole in the sea quark sector {R'jj^i) vanishes, and the re- 
maining ones satisfy a larger number of relations than for 
rfsoa = 3. The remaining nontrivial residues for dsca = 2 
may be expressed as 



R]k 




Rf 




Rl 


= R\jj + 



R, 



6jj 



Til 



(9) 



It is also noteworthy that the above residue notation is 
highly redundant because of the trivial relation between 
X47X6 and Xv This fact has been exploited in the sim- 
plification of the end results. 

For dsoa = 3, the naturally generated residues are 



R 



■jkl 

Rf 
Ri 



i456jkl 1 



= R 



Til _ TDi TDl 

ixq^^-j xi.^j^j^ n^j^^ . 



(10) 



Even in this case there are many relations between the 
various residues. Some can be found in Ref. 4] but many 
more exist. Finally, it should be noted that in the limit 
where the sea quark masses become equal to the valence 
quark masses, the propagator residues of PQxPT reduce 
so that the ir^ and ry meson propagators of unquenched 
xPT are recovered. 



D. Integral Notation 

The expression for the decay constant of a charged 
pseudoscalar meson depends, in NNLO PQxPT, on a 
number of one-loop and two-loop integrals. The finite 
parts of the chiral logarithms A, B and C are 



Mx) 
c(x,x,x;0) 



7ri6Xlog(x/M^), 

Xi log(Xi/M^) - Xj log(Xj7M^) 



Xt - Xj 



-7ri6/(2x), 



(11) 



where the subtraction scale dependence has been moved 
into the loop integrals. We also define ttiq = l/(167r^). 
Note that in the limit Xi — Xji the integral B reduces to 

^(X,X;0) = -7ri6(l + log(x/M2)). (12) 
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The following combinations of finite one-loop integrals 
are also introduced, as they are naturally generated by 
the procedure of dimensional regularization: 

A[x\e) = A(x)V(27ri6x) 

+ 7ri6x(^Vl2 + l/2), 

B(x,x;0,e) = A(x)B(x,X;0)/(7ri6x) 

- A(x)V(27ri6x') 

+ vrie (^Vl2 + 1/2). (13) 

Additionally, the upper middle and upper left diagrams 
in Fig. n generate integrals of the form 



(14) 



which are symmetric under the interchange of Xi s-nd 
Xj ■ The above integrals have been introduced in order to 
make the symmetries in the end results more explicit. 

The finite two- loop integrals H^,H[,H2i that are 
generated by the top right diagram of Fig. may be 
evaluated using the methods of Ref . • Note that the 
corresponding primed integrals , Hi , H21 indicate 
differentiation with respect to . The notation used for 
the H integrals in this paper is similar to that of Ref. ^3 j 
except that an extra integer argument now indicates the 
propagator structure, such that e.g. 



H''iXi,Xj,Xk;p'^) H'^{n, Xi,Xj,Xk]P^)- 



(15) 



The case of n = 1 indicates that the integral consists of 
single propagators only, as in Ref. [l^l , whereas n = 2 
indicates that the first propagator appears squared and 
n — i that the second propagator appears squared. The 
cases with two double propagators that can appear in 
the calculations are n = 5, for which the first and second 
propagators appear squared, and n = 7 for which the sec- 
ond and third propagators are squared. Explicit expres- 
sions for n — 1 can be found in Ref. 13], and the other 
cases may be obtained by differentiation with respect to 
the masses of those expressions. 



Xbdn 



-ace 
Abdn 



3 ^ 

i = 4,5,6 



3 ^ 

i = 4,5,6 



DO pc n 
^bi ^di Xi 5 



Dti pc pe n 
^bi ^di ^ii Xi 



pi I pi rj pz 



(17) 



For dsca = 3 these become 



Xbc 



Xbcefn 



-adk 
Xbceflr, 



TfV 

^ijkl 



^bii Xi ) 

i = 4,5,6 

^bci Xi 5 

i = 4,5,6 

X! ^bciRifiXi, 
i = 4,5,6 



\ ^ pa pd pfe n 

o "-bci ^efi ^lii Xi ) 



1 

3 

i = 4,5,6 

pi 1 pi n pi 

^]kk + ^ ^^jfei- 



(18) 



It should be noted that also in the case of dgea = 2, the 
sums run over all three sea quark flavors. Furthermore, 
only some of the above x functions appear in the results, 
as many of them can be reexpressed in terms of other R 
or functions. 

The actual simplification of the end results is at first 
glance a formidable task, since the expression for Fa as 
calculated from the diagrams in Fig. ^ contains several 
thousand terms. Especially for dvai — 2, the parts pro- 
portional to yl(xj);£) and ^(xi;e) alone are several hun- 
dred terms long, while in simplified form they typically 
contain only ~ 5 terms. However, these simplifications 
are not easily apparent and are consequently best accom- 
plished by the employment of suitable software, such as 
Maple or Mathematica. Even so, as the expressions can 
seldom be factored into a single term, considerable trial 
and error is usually required. 



E. Simplification of End Results 

A number of combinations of quark masses and propa- 
gator residues are naturally generated in the calculations, 
and have consequently been given special notations. The 
most common of these is 

i = 4,5,6 

of which xi is equal to the average sea quark mass x 
defined in Ref. Other combinations that appear in 
the calculations consist of products of quark masses and 
propagator residues. For dgca = 2, these include 

Xbn ^ ^ ^ ^bi Xi t 

i = 4,5,6 



F. Summation Conventions 

The expressions for the decay constants are symmetric, 
within the sea and valence quark sectors, under the in- 
terchange of quark masses. Consequently, the end results 
can be conveniently compactified by the introduction of 
summation conventions which exploit these symmetries. 
For example, the notation for the sea quark sector may 
be considerably compactified by the introduction of the 
summation indices s and t, which may appear for the 
quark masses x? and among the indices of the functions 
R and x- 

These sea-quark summation indices should be inter- 
preted as follows: If an index t is present once or several 
times, then there will always be an occurrence of the s 
index as well, and the term is then to be summed over 
all pairs of different sea quark indices. If the index s is 
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present but t is not, then the entire term is to be summed 
over all sea quark indices. Simple examples of terms are 



III. THE PSEUDOSCALAR MESON DECAY 
CONSTANTS TO 



Xst 



1 = 4,5,6 

i,j = 4,5,6 
j>i 



(19) 



For terms consisting of products of several factors, the 
summation sign should always be inserted at the begin- 
ning of each relevant term, such that e.g. 



rji Xi ■ 



(20) 



The decay constants Fa of the pseudoscalar mesons are 
obtained from the definition 



{0\A^mMp))=^V2p^^Fa, 



(23) 



1 = 4,5,6 



in terms of the axial current operator A^. The diagrams 
that contribute to that operator at 0{p^), or NNLO, are 
shown in Fig. ^ Diagrams at 0{p'^) and 0{p'^) also con- 
tribute to Eq. idSl via the renormalization of the pseu- 
doscalar meson wave function (j)aip)- The results for 
charged pseudoscalar mesons so obtained depend on the 
0{p^) low-energy constants if[g through i^as. 



Thus any contribution which is written in terms of the 
indices s and t explicitly fulfills the required symmetry 
properties in the sea quark sector. Note that also for 
dsoa = 2, the summation is over all three sea quark fla- 
vors, although the index t has not been implemented in 
that case. 

Further compactification of the results is possible, as 
the valence quark sector is symmetric under the exchange 
of the valence quark indices 1 and 3. For this purpose, 
the summation indices p and q have been introduced. 
These are in the present paper only needed for the case 
of dvai = 2 and then always occur for the squared valence 
quark masses xi f^^d xs- If the index q is present, there 
will always be an index p and the resulting sum is over 
the pairs {p,q) ~ (1,3) and {p,q) = (3,1). If only p is 
present, the sum is over the indices 1 and 3. An example 
is 





FIG. 1: Feynman diagrams at 0{p^) or two-loop for the 
matrix element of the axial current operator. Filled circles 
denote vertices of the £2 Lagrangian, whereas open squares 
and shaded diamonds denote vertices of the £4 and La La- 
grangians, respectively. 



MXp)R^gnXp = MXl)Rlr,Xl 



(21) 



Any contribution written in terms of the (p, q) notation 
is thus symmetric under the interchange of the valence 
quark masses xi and xs- 

For dson = 3, there exists an additional symmetry, i.e. 
the results are symmetric under the interchange of the 
lowest order neutral meson masses, Xtt and Xri, in the sea- 
quark sector. This is exploited by the summation indices 
m and n. If the index m is present, there will always be 
an index n and the corresponding sum is over the pairs 
{m,n) = {TTjrj) and (rn,n) = (?7,7r). An example is 



The results are expressed in terms of the valence inputs 
Xi, X3 and the sea inputs X4j X5j Xe, which are defined in 
terms of the quark masses through Xi = '^Bomi, and 
the quantity Xij — iXi + Xj)/2, which corresponds to 
the lowest order charged meson mass. Other parameters 
include the decay constant in the chiral limit (Fq), the 
quark condensate in the chiral limit, via (qq) = —BqF^, 
and the LEC:s of 0{p^) and 0{p^), i.e. the L\ and the 
[Hill, respectively 

The decay constants of the pseudoscalar mesons are 
given in the form 



A{Xm)KiiXm = A{x^)R!l^^Xn + ^ A- (22) 

The above summation conventions thus provide a means 
of eliminating a large number of terms in the expression 
for Fa, which are of similar form. As an added bonus, it 
gives a possibility to conveniently check that the NNLO 
loop results fulfill the required symmetry relations in the 
sea and valence sectors. 



phys 



^(4)vi 
-^0 



c(6)vs 



-I- 



r(6)vs 
loops 



0{p^ 



(24) 

where the 0{p'^) and 0{p^) contributions have been sepa- 
rated. The NNLO contribution S'^^^ has been further split 
into the contributions from the chiral loops and from the 
0{p^) counterterms. The superscripts (v) and (s) indi- 
cate the values of dvai and c?soa, respectively. 
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A. Results for dvai = 1 

The NLO result for dvai = 1 is fairly short, and will 
thus only be given for dgea — 3. The results for dgcn — 
1,2 may readily be derived from that expression. The 
combined NLO result (loops and counterterms) for dsca — 
3, which is in agreement with Ref. 0, is 

JW13 ^ \2L\xi + 4L^xi + l/2A(xi.). (25) 

At NNLO, the chiral loops form, by far, the largest 
contribution to the decay constant. As a straightforward 
derivation of the results for dgea = 1:2 from the dgoa = 3 
case is tedious and complicated, the expressions for the 
different cases are given separately below. As expected, 
the infinities in all expressions for the decay constant have 
canceled. The appearance of 'unphysical' B logarithms 



in the results is, in part, due to the partial quenching, and 
to the fact that the results have been expressed in terms 
of the lowest order masses rather than the full physical 
ones. Consequently, not all of them correspond to the 
quenched chiral logarithms which are ill-behaved in the 
chiral limit. 

The contribution from the Oijf'^ counterterms to the 
decay constant at NNLO is, for dgoa = 3, 

^ - 8X[9X? + 24^2^0X1X1 + 24i^2^iX2 
+ lIKl^xl + 8X2^3X?, (26) 

from which the results for c?sca = 1,2 may be readily 
inferred. The NNLO contributions from the chiral loops 
for dsca = 1,2,3 are, respectively. 



Cops = ^16 [-13/3x1X4 + 1/2X?- 3/2x1] - 2 7ri6L^x? - ^i^^l [X1+8X4] 

+ ^leiS [-17/6x1X4 + 5/4X?- 3/4x1] + vr^g [-1/2 xiX4 + 1/64 x? - 73/128 x^] 

- ^%L\LlxiX^ - 72Lf xl - 8i^^x? + 4A(xi)£S [xi + i?f] - 4A(xi)L^Xi - 10A(xi)£^Xi 
+ 4A(xi)LS [Xi+^f] - 4/3A(xi)L^Xi - 1/2 A(xi)B(xi4, Xw; 0) Xi4 + 1/8 A(xi; e) ttis X4 

- ^(Xi4)vri6 [3/4x14 + 9/8x4] - 12 l(xi4) Xi4 - 30 ^(xu) Xi4 - 18A(xi4)i4X4 

+ 6l(xi4)iSxi + 9/4 A(xi4;e) ^16X4 - MA{x^)L\x^ - 16A(x4)iSx4 + 32A(x4)i^X4 

+ A(x4;e)^i6X4 + 4B(xi,Xi;0)L5i?fxi + 4B(xi,Xi;0)L5^fxi - 4/3 5(xi,Xi;0)LEi?fxi 

- l/8B(xi,Xi;0,e)^i6i?fxi - 36 5(xi4,Xi4;0)L^Xi4X4 - 12 B(xi4, Xw; 0) x?4 

+ 72B(xi4,Xi4;0)Lgxi4X4 + 24B(xi4,Xi4;0)igx?4 - 1/8 i/^(l, xi, Xi4, Xw; Xi) [xi-^f] 

- ■H"^(i,xi4,xi4,X4;xi)x4 + i/8ff^(2,xi,xi4,xi4;xi)-RiXi + 5/i8iJ^ (i,xi,xi,xi;xi)xi 
+ -ff^ (i,xi,xi4,xi4;xi) [1/8x1X4 - 1/2 xi] + -ff^ (i,xi4,xi4,X4;xi)xiX4 

+ 2/9i7^'(2,xi,Xi,Xi;Xi)^?X? + 3/8iI^'(2,xi,Xi4,Xi4;xi)i??X? + i/9i7^'(5,xi,xi,xi;xi)(^?)'x? 

- 2fff (3,xi4,Xi,Xi4;Xi)^iXi + 3/8ff|i (l,xi,Xi4,Xi4;Xi)Xi + 3i?fi'(l,X4,Xi4,Xi4;Xi)Xi 

- 3/8i?fi'(2 , Xi,xi4,xi4; Xi)i??Xi, (27) 

Cops - ^16 [4/9 Xr,X4 + 1/2 x\ - 13/3 XiXi - 35/18 X2] - 2 L\ x\ 

- [11/3 Xr,X4 + Xi + 13/3x2] +^i6iS [4/9 x^X4 + 5/4 X?- 17/6x1X1 -43/36x2] 

+ [- 15/64 x„X4 + 1/64 X?- 1/2x1X1 -43/128x2] - 48i^L^xiXi - 72if x? - 8if Xi 

+ 4l(x,)LS^?iX, - 8A(x,)L^Xr, - 2A(x,)L^X, + 4 A(x,) ^?i X, +4A(x,)i4X, 

- 4/3A(x^)L^i??iXi - l/6A(x^)B(xi.,Xi.;0)i?Lxi. + 1/8 i(x„; e) ^le i?J Xr, 

+ 4l(xi)L5 [i??Xi+^f] -4A(xi)£^Xi - 10A(xi)LSxi +4A(xi)LS [i?J Xi + ^f] 

- 4/3l(xi)LSi?Jxi - l/6A(xi)B(xi.,Xi.;0)i?,\xi. + A(xi;£)^i6 [l/4xi - l/8i?Jxi - l/8i?f] 

- 24A(x4)iIX4 - 6A(x4)i2X4 + 12A(x4)i4X4 + 3/8 l(x4;e) ^16X4 - 32A(x46)i^X46 

- 8l(x46)i2X46 + 16A(x46)i4X46 + 1/2 A(x46; e) TTie X46 - A(xi.) TTig [1/4 Xls + 3/8 Xl] 

- 4l(xi.)ioXi. - 10l(xi.)iSxi. - 6l(xi.)i^Xi + 2A(xis)i5Xi + 3/8 l(xi.; e) ^le [Xs + Xi] 
+ 4B(xi,Xi;0)iS^fxi + 4B(xi,Xi;0)£5i?^Xi - 4/3 5(xi,Xi;0)i^i?fxi 

- l/8B(xi,Xi;0,£)^i6i??Xi - 12B(xi.,Xis;0)L^XiXi. - 4B(xi.,Xi.;0)L5xL 

+ 24B(xi.,Xi.;0)i6XiXis + 8B(xi.,Xi.;0)L8xL + l/24i/^(l,x^,Xis,Xis;Xi)^:;i.Xr, 



+ H^{l,Xi-Xu,Xu;Xi) [- 1/12 ii^^ XI + 1/24 XI + 1/24 iif] - 3/8 F^(l, Xi4, Xi4, X4; Xi) X4 

- l/2i/^(l,xi4,Xi6,X46;Xi)X46 + l/24ff^(2,xi,Xi^,Xu;Xi)^fxi 

+ i/9F^'(i,x^,Xn,xi;xi)Ki)'x? + 2/9ff^'(i,x^,xi,xi;xi)i??ii?Jx? 

- i/^'(l,X^,Xi.,Xi,,;Xi) [l/6i??iX? + l/24i?:;i,XnXi] + i^^'(l,Xi,Xi,Xi;Xi) [i/6x? 
+ 1/9 (R'.fxl] + i/^'(l,xi,Xi.,Xis;Xi) [-1/4 i?,\xi + 1/8 i?Jx?- 1/24 J?^Xi] 

+ 3/8ff^ (l,Xl4,Xl4,X4;Xl)XlX4 + l/2i?^'(l,Xl4,Xl6,X46;Xl)XlX46 

+ 2/9F^'(2,xi,Xr„Xi;Xi)i??i^'?X? + 2/9ff^'(2,xi,Xi,Xi;Xi)^5^fx? 

+ l/8ff^'(2,xi,Xi.,Xi.;Xi)^'fx? + i/9-ff^'(5,xi,xi,xi;xi)(^f)'x? 

- i/3Jff' (i,xr„xis,xi.;xi)^'L^L^x? - 2/3Jff' (i,xis,xi.,xi;xiX^Lr,x? 

- 2/3i/f (3,xis,xi,xis;xi)-Ri Xi - i/8i?fi(i,Xr,,xi,9.xis;xi)^r"isXi 

+ if|;'(l,Xi,Xi.,Xi.;Xi) [l/4i?,i^x?-l/8i??X?] + 9/8i/2i'(l.X4,Xi4,Xi4;Xi)x? 
+ 3/2ijf;(l,X46,Xi4,Xi6;Xi)x? - l/8i/2''i'(2,Xi,Xi«,Xi.;Xi)i?fx?, 

= Trie [4/9 X,.Xr, + 1/2 X? - 13/3 XiXi - 35/18 X2] - 2 L{ xl 

- ^16 [11/3 x.X^ + X? + 13/3x2] +^16^5 [4/9 x.Xr, + 5/4 X?- 17/6x1X1 -43/36x2] 

+ [-15/64 x.Xr, + 1/64 X?- 1/2x1X1 -43/128x2] - 48i^L^xiXi - 72 x? " 8Ll''xl 
+ 4 Aixm) LI i?™ 1 Xm + 8 A(x™) L[ X™o X™ + 2 A(x„0 ^3 x"o X™ + 4 ^(x™) ^3 ^^ii Xm 

- AA{xm)LlXni - 4/3l(x™)L^K'\iXi " 1/6 A(x„)S(xi„ Xi., 0)K\, xi. 

+ A(x™;e)7ri6 [-7/24 x™o Xm + 1/6 x"i - 1/8 ii™iiXm] +4A(xi)i.S [R\xi+Ri] 

- AA{xi)Llxi - 10A(xi)i.^Xi +AA{xi)Ll [Rlxi+Ri] - A/Z A{xi) LlRlxi 

- l/6i(xi)B(xi«,Xi«,0)iiUxi. +i(xi;£)7ri6 [l/4xi - l/8i?? Xi - l/8i?f] 

- ^(xi.)vri6 [1/4X1S + 3/8x1] - 4l(xi.)i:5xi. " 10 A(xi.) Xi. - 6A(xi.)LIxi 
+ 2A(xi.)L;;xi + 3/8A(xi.;e)7ri6 [xs + Xi] + Mxs) LI [-8 x. + 8/3 i?^ Xs] 

- A{xs)Ll [2xs-2/3i?^x.] + Mxs)Ll [4xs - 4/3i?:x.] + Mxs\e)^^Q [1/8 x. - 1/24 i?^ X.] 

- 1& A{xst) Llxst - AA{xst)L^2Xst + 8A{xst)Llxst + l/4A(x.t;e)7ri6Xst 
+ 4B(xi,Xi;0)LS^fxi + 4B(xi,Xi;0)i5^«?Xi - 4/3B(xi,Xi;0)LSi?fxi 

- l/8B(xi,Xi;0,e)7ri6i??Xi - 12 B(xi., Xi^; 0) Xi Xi« " 4 B(xi., Xi«; 0) xL 

+ 24B(xi.,Xi.;0)L;5xiXis + 8 B(xi., Xi.; 0) xL + l/24i?^(l,Xm,Xi.,Xi.;Xi)^™„isXm 

+ //^(l,Xi,Xi.,Xis;Xi) [- 1/12 i?i,^xi + 1/24 i?Jxi + 1/24 i?f] + i?^(l, Xi., Xis, x«; Xi) [-1/8X. 

+ l/24i?,= xs] - l/4i/^(l,Xi,s.Xu,X.t;Xi)x.t + l/24i/^(2,xi,Xi.,Xis;Xi)^^fxi 

+ l/9iJ^'(l,Xr„,Xm,Xi;Xi)«'ii)'x? + 2/9//^'(l,Xm,Xi,Xi;Xi)Kn^?X? 

- if^'(l,Xm,Xi«,Xi.;Xi) [l/6KnX? + l/24i?:;,„i,XmXi] 

+ 2/9H^'{l,x.,Xr,^Xi;Xi)R;uKnXl + if^'(l,Xi,Xi,Xi;Xi) [i/6x? + i/9(J??)'x?] 

+ if^'(l,Xi,Xia,Xia;Xi) [- 1/4 i?L^x? + 1/8 i^Jxf- 1/24 iifxi] + -^""(1, Xi«, Xi«, X«; Xi) [1/8X1X. 

- l/24i?^xiX«] + l/4ff^'(l,Xi«,Xu,Xat;Xi)XiX«t + 2/9ff^'(2,xi,Xm,Xi;Xi)J«niififx? 

+ 2/9H'"{2,xi,xi,xi;xi)RiRtxi + {2,xi,xis,xu;xi) Rfxi 

+ 1/9H^' {^,xi,Xi,Xi\Xi) {Rifxl - l/3H['{l,Xm,Xu,Xu;Xi)RnuRlsmXi 

- 2/SH['{l,xu,Xis,Xi;Xi)RnisRumXl - 2/SH['{3,xu,XuXu;Xi)Rfxl 

- ^/»Hi|{l,Xm,Xls,Xls■,Xl)Rl^nuX^ + H[^ii,xi,Xis;xis;xi) [i/4i?Lr,x?-i/8i??x?] 
+ i?2''i'(l,X«,Xia,Xi«;Xi) [3/8x?-l/8i?^X?] + 3/4iff;(i,x.t,xi«,xit;xi)x? 

- l/8J?2''i'(2,Xi,Xi«,Xi«;Xi)^iX?- 
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B. Results for dvai ~ 2 

In general, the expressions for dvai ~ 2 are longer than 
those for dvai = 1, as the number of independent integrals 
that can appear is significantly larger. The full result for 
dsca = 3 at NNLO is very large, and has not yet been 
worked out. Consequently, the largest number of nonde- 
generate sea quarks considered for dvai = 2 is c?sca = 2. 
The NLO result for rfvai = 2 and dsca — 2, which agrees 
with the result of Ref. p\ is 

(5(4)22 ^ ^2L^xi + 4L^Xi3 + Aixp) [1/6 

- 1/12 i?^] + l/4A{xps) - l/12A(x,)i?;i3 



- l/12B{xp,Xp,0)Rt- (30) 

At NNLO, the contribution to the decay constant from 
the 0{p^) countcrtcrms for rfvai = 2 is similar to that for 
t^vai = 1- For dsca = 2, that contribution is 

' = AKl.xl + 24^2^0X1X13 + 24/^2^1X2 
+ 72 K-,^xl + 8K;,xiX3- (31) 

The NNLO contributions to the decay constant from 
the chiral loops are, for dsea = 1 and dsea = 2, respec- 
tively. 



Cops = 7^16 [-1/2 X1X3- 13/3 xi3X4 + X?3- 3/2x4] - 2n^,Llxl3 ~ ^16^ [xls + ^xl] 

+ TTieLJ [-7/12 X1X3 -17/6x13X4 + 11/6 x?3- 3/4x4] + [-59/384 xiX3 - 1/2 Xi3X4 + 65/384 xL 

- 73/128 x1] - 48L^L^xi3X4 - 72Lf x' - 8if x?3 + A{Xp)t^w [5/96 Xp + 1/32 Xg + 1/48 X4 

- 1/12 i?^xp- 1/8 i??X4] +AiXp)L^o [2/3 xp + 8/3 i?Pxp + 2/3 i?^] + A{xp) [5/3 Xp + 2/3 i?^ Xp 
+ 5/3 Rf] +A{xp)Ll [x4-2i?Px4] +A{xp)Ll [I/6 Xp - 1/6 X9 - 1/3 X4 - 2/3 i?^ Xp] 

+ Aixpf [19/288- 1/72 i?Pi?^] - A(xp)A(xp4) [I/I6 + 1/12 i?^] - A{xp)A{x,4) [3/16 + 1/12 i?^] 
+ l/8i(xp)A(xi3) + i(xp)B(Xp,Xp;0) [1 1/36 Xp- 1/18 i?Pxp- 1/72 + 1/144 i?3 

+ i(xp)B(x„X9;0) [-l/72i?Pi?^ + 1/144 i?^] - l/4i(xp)S(xp4,Xp4;0)xp4 

- l/18A(xp)B(xi,X3;0)i?^Xp + l/18l(xp)C(Xp,Xp,Xp;0)i?^Xp + l/8A(xp;£)^i6 [-Xp + i??X4] 

+ ^(Xp4) 71-16 [3/16 Xp4 - 9/16 Xg4 - 9/16x4] - 6^(Xp4)^SXp4 - 15^(Xp4)i5Xp4 - 9^(Xp4)i4X4 

+ 3i(xp4)i5Xi3 - 9/32A(xp4)' + i(xp4)^(Xp,Xp;0) [1/8 Xp - 5/8 Xp4] - 1/16 A(xp4)S(x„ X,; 0) i?^ 

+ l/6A(xp4)5(xi,X3;0)x4 + l/3^(xp4)S(xi,X3;0,fc) + 9/8 A(xp4;£)^i6X4 + A(xi)A(x3) [-1/144 

+ 1/36 RlRl] - 4i(xi3)iJxi3 - 10i(xi3)^^Xi3 + l/8A(xi3)' - 1/2 i(xi3)S(xi, X3; 0, fc) 

+ 1/4 A(xi3;e) ^16X13 + 9/16i(xi4)^(x34) - 64i(x4)^^X4 - 16A(x4)i2X4 + 32A(x4)i4X4 

+ 2/9^(x4)5(Xp,Xp;0)x4 - 4/9^(x4)S(xi,X3;0)x4 + A(x4;e)^i6X4 + S(xp,Xp;0)^i6 [l/96i?;^Xp 

+ 1/32 i?^X9 + 1/16 i?^X4] + 2/3B(xp,Xp;0)L5i?^Xp + 5/3S(xp,Xp;0)2.S^^Xp 

+ B(Xp,Xp;0)i^ [2xpX4-4i?PxpX4 + 3i?^X4] + S(xp, Xp; 0) [-2/3 XpX4 + 4/3 Xp - 4/3 i?^ XpXi3 

- l/3i?^xi3] + B(Xp,Xp;0)L^ [4x4-8i?^XpX4] + 4 B(xp, Xp; 0) L^^p)' + B(Xp, Xp; 0) [4/3x1 

- 8/3 i?«Xp] +^(Xp,Xp;0)2 [1/18 i?«i?^xp + 1/288 (i?^)2] - 1/18 B(xp, Xp; 0)B(xi, Xs; 0) ii^i?^ Xp 
+ l/18 5(xp,Xp;0)C'(xp J Xp I Xp ; 0)(i?^)'Xp - 1/8 B(xp,Xp;0,e) Trie i?^Xp4 

- 18i?(Xp4,Xp4;0)L4Xp4X4 - 6i?(Xp4,Xp4;0)L5Xp4 + 36i?(Xp4,Xp4;0)L6Xp4X4 

+ i2 5(xp4,Xp4;0)L^Xp4 + 1/144 s(xi, XI ;0)B(x3,X3;0)i?fi?3^ - 8S(xi,x3;0)L;i?^i?^ 

- 8/3 5(xi,X3;0)i^i?fi?3^ + 4C(xp,Xp,Xp;0)i^i?^XpX4 + 4/3C'(xp,Xp,Xp;0)i^i?^Xp 

- 8C(xp I Xp I Xp ;0)igi?pXpX4 - 8/3 C(xp I Xp 7 Xp ; I Xp? Xpi Xi3j X13) [1/8 Xp 

- 5/72 x13 + 1/36 i?^i?«xi3] + i^''(l,Xp,Xi4,X34;Xi3) [-1/16 Xp + 1/16 Xg - 1/8 i?? X4] 

+ i/^(l,Xl,Xl3,X3;Xl3) [1/72x13- 1/18 i?Ji??Xl3] - l/8i/^(l,Xl3,Xl3,Xl3;Xl3)Xl3 

- ^f^(i,xi4,X34,X4;xi3)x4 + i^^(2,Xp,Xp,Xl3;xl3) [i/36 i?^i?p xi3 - i/72i?pXi3] 

+ i/^(2,Xp,Xi3,X9;Xi3) [l/36i?«i?^xi3-l/72i?^Xi3] + 1/8 i?^(2, Xp, Xi4, X34; X13) i?p Xp4 

- 1/144 i/^(5,Xp,Xp,Xi3;Xi3)(i?^)'xi3 - l/72i/^(5,xi,X3,Xi3;Xi3)i?^i?^Xi3 



9 



+ H^'{l,Xp,Xp,Xi3;Xi3) [-l/8XpXi3 - 5/24x?3 + l/12i??X?3 - 1/36 (i?^' xL] 

+ H^' (1, Xp, X14, X34; Xis) [3/16 XpXi3 - 1/16 Xi3X4 - 1/8 RpR^ Xi3 + 3/8 i?« x^s] 

+ i?^'(i,xi,Xi3,X3;xi3) [-1/72x13 + 1/QRlRlxh] + i/8ff^'(i,xi3,xi3,xi3;xi3)x?3 

+ i?^'(l,Xi4,X34,X4;Xi3)Xi3X4 + ^f^' (2, Xp, Xp, Xi3; X13) [-l/36i?Pi?;^x?3 + 5/72i?;!x?3] 

+ //^'(2.Xp,Xi3,X3;Xi3) [-l/36i?«i?^x?3 + l/72i?^X?3] - 1/8 iJ^'(2, Xp, Xi4, X34; Xw) i?^ Xp4Xi3 

+ 5/144 i?^'(5,Xp,Xp,Xi3;Xi3)(i?^)'x?3 + l/72ff^'(5,xi,X3,Xi3;Xi3)i?^i?^X?3 

+ iff (l,Xp,Xp,Xi3;Xi3) [10/9 x?3- 2/9 iJ^iJ'xL] - ^fr(l,Xp,Xi4,X34;Xi3)i?^X?3 

- ^fr(i,Xp4,x,4,Xp;xi3)x?3 + ^^r(i,xi3,xi,X3;xi3) [i/9xL - 2/9i?^i??x?3] 

+ ^fr(3,Xi3,Xp,Xp;Xi3) [l/9i?Pi?^x?3-l/9i?^X?3] + l/9iJr{3,Xi3,Xp,Xg;Xi3)i?^i?^X?3 

- l/18i/r(7,Xi3,Xp,Xp;Xi3)(i?p)'x?3 - 3/8i72^i'(l,Xp,Xp,Xi3;Xi3)x?3 
+ 3/8 Hi, {1 

, Xp, Xi4, X34; Xi3)i?^?X?3 + 3/8 Hi {1 

7 Xp47 Xpi Xq4: ;xi3)^gX?3 

- 3/8i?2';'(l,Xp4,X9,X94;Xi3)i?^X?3 + ^^l^i' (1 , Xi3 , Xp, Xp i Xi3 ) [5/24x?3-l/12i?^i?^X?3] 

+ i^2''l'(l,Xl3,Xl,X3;Xl3) [-l/24x?3 + l/6i?3^?X?3] + 3/8ifil'(l,Xl3,Xl3,Xl3;Xl3)X?3 

+ 3 7?|; (l,X4,Xl4,X34;Xl3)Xl3 - 3/8i?|; (3,Xp4,Xp,X<?4;Xl3)-RpXl3 

+ i?|'i'(3,Xi3,Xp,Xp;Xi3) [-l/12i?^i?^x?3 + l/24i?^X?3] + ^^i^i' (3, Xi3, Xp, X,; X13) [-l/12i?«i?^x?3 

+ l/24i?^xi3] + l/48i?i'i'(7,Xi3,Xp,Xp;Xi3)(i?^)'x?3 + 1/24 i?i'i'(7, Xi3, Xi, X3; X13) X?3, (32) 

= ^16 [4/9 X^X4 - 1/2 X1X3 + X?3 - 13/3 X1X13 - 35/18 X2] - 2 w,e L\ x?3 

- TTiei^ [11/3 x^X4 + Xi3 + 13/3 X2] + Trie [4/9 x„X4 - 7/12 xiXs + 11/6 X13 - 17/6 X1X13 - 43/36x2] 
+ TT^e [- 15/64 x^X4 - 59/384x1X3 + 65/384 x?3- 1/2X1X13 -43/128x2] - 48L^L^xiXi3 - 72if x? 

- 8Lfx?3 +^(Xp)'ri6 [-1/24 Xp + 1/48 xi- 1/8 Xii?^^ + 1/16 Xii?p- 1/48 i?^^Xp- 1/16 i?^,Xg 
+ l/48i?^pXn + l/16i?^Xi3] + A{Xv)Ll [8/3 ii^^ Xp + 2/3 iJ^ Xp + 2/3 i?^] + A{xp) [2/3 i?^^Xp 
+ 5/3 Rl Xp + 5/3 R^] + A{xp) LI [-2 xiX^^o -2xiR% + 3 xi R'p] + A(xp) LI [-2/3 x^i - Xp 
+ 1/3 RP^ X, + 1/2 R; Xp " 1/6 R; X,] + ^(Xp)' [I/I6 + 1/72 (RP^f - 1/72 RP^R; + 1/288 [R^f] 

+ Aixp)A{xps) [-1/36 RP,^- 5/72 RP^ + 7 /UAR;] - A{xp)A{x,s) [l/36 i?^^ + 1/24 i??,, + 1/48 E^] 

+ A{xp)A{xv) [-1/72RP^R'^,, + 1/UAR;r;,,] + l/8A(xp)i(xi3) + 1/12 A(xp)i(x46) i?^p 

+ A{xp)Bixp, Xp; 0) [1/4 Xp - 1/18 RP^R; Xp " 1/72 RP^R"^ + 1/18 [R^f Xp + 1/144 R^R^] 

+ A{xp)B{xp,Xv;0) [1/18 i?^j,i?^Xp- 1/18 i??3^^Xp] + A(xp)5(x„Xg;0) [-1/72 i?^X + 1/144 i?X] 

- i/i2A(xp)5(xp.,Xp.;0)i??^Xp. - i/i8A(xp)5(xi,x3;0)i?«^i?;xp 

+ 1/18 A{xp)C{xp, Xp, Xp-, 0) R;Rt Xp + A{Xp; £) Trie [I/8 XiR^r, ' l/16Xi R^ - 1/16 R; Xp - 1/16 R^] 

+ ^(Xps)7rie [1/16 Xps - 3/16 Xqs - 3/16x1] - 2A{xps)Llxps - 5A{xps)Llxps - 3^(xps)i4Xi 

+ A{xps) LI X13 + i(xp,s)^(x,) [7/144 R^^ - 5/72 R;, - 1/48 R^ + 5/72 R^ - 1/36 R^^] 

+ A{xps)Bixp,Xp;0) [l/24i?P^Xp- 5/24 i?P^Xp.] + A(xp.)B(Xp, Xr,; 0) [-1/18 i?^,i?^j,^Xp 

- l/9i?^,i?^^p^xp.] - l/48A(xp.)B(x„X9;0)i?^ + l/18A(xp.)S(xi,X3;0)i?,%x. 

+ l/9i(xps)^(xi,X3;0,/c)i?,% + 3/16A(xp.;£)^i6 [x. + Xi] - l/8^(Xp4)' - 1/8 i(xp4M(Xp6) 

+ l/8^(xp4)^(x96) - l/32A(xp6)' +^(x^)^i6 [1/16 xii?:;i3- 1/48 i?:;i3X^ + 1/16 i?;;i3Xi3] 

+ A{x^)U, [4R'I,x, + 2/3R;,3Xv] - 8A{x^)Llxr, - 2A{xn)Llxn + Mxn) [AR\^Xv + h/i R^^^Xv] 

+ Aix^)Ll [4x, + Xii?;;i3] - Mxv)Ll [l/6i?^pX, + ^?3Xi3 + l/6i?;i3X^] + 1/288 l(x^)2 (i?;i3)2 

+ l/12A(x^)A(x4e)ii;i3 + ^(Xr,)B(Xp,Xp;0) [-1/36 x^^? - 1/18 iJ^^iJ^^Xp + 1/18 fl;ipV?;5Xp 

+ l/144i?Xi3] +^(Xr,)S(Xp,X.,;0) [-l/18x^^K 1/18X^^? + 1/18 Wi?^p,Xp] 
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- 1/12 A(xr,)B(xps,Xps;Q)RlsXps - A(x^)B(x,„x„;0) [l/216ii;i3X4 + l/27i?,;;i3X6] 

- l/18A(x^)B(xi,X3;0)i?i^i?^^X^ + l/18A(x^)(7(xp,Xp,Xp;0)i?^Xxp + ^^Xn^ 

e)7ri6 [1/8 x»j 

- 1/I6xii?;i3 - l/8i?l'3X^ - l/16i?;i3X^] + A(xi)A(x3) [-1/72 Rl^Rl + l/mRl^R% + 1/144 

- 4^(xi3)i'iXi3 - 10^(xi3)i^Xi3 + l/8A(xi3)' - l/2A(xi3)S(xi,X3;0,fc) 

+ l/4i(xi3;£)7ri6Xl3 + l/4^(Xl4)v4(X34) + l/16^(xi6)^(_X36) - 24^(x4)i'iX4 - &A{xA)Llxi 

+ 12A{xi)Llxi + l/l2A{xi)B{xp,Xp\^){R%fXi + l/6^(x4)B(Xp,Xr,;0) [R^R^^Xi - Rl^R'^AXi] 

- l/24A(x4)-B(x^,x„;0)i?;i3X4 - l/6A(x4)S(xi,X3;0)iii^i?!^X4 + 3/8A(x4;£)7ri6X4 

- 32A(x46)iiX46 - 8l(X46)i2X46 + 16 1(X46) ^4 X46 + A{Xie)B{Xp,Xp\^) [l/9 X46 + 1/12 ii^p Xp 

+ l/36i?^pX4+ l/9i?^4X6] + A(x46)S(xp,Xr,;0) [-1/18 ii^^ X4 - 1/9 ii^4 X6 + 1/9 i?^4 X6 + 1/18 i??3 X4] 

- i/6i(x46)B(xp,x.,;0,fc) [R^p-RU + i/9A(x46)s(x„x^;0)i?;i3X46 - ^(x46)s(xi,X3;0) [2/9x46 

+ l/9i?;;4X6+ l/18i??3X4] - l/6i(x46)S(xi,X3;0,fc)i?l'3 + 1/2 i(x46 i e) TTig X46 

+ S(xp,Xp;0)7ri6 [1/16 xii?^ + 1/96 i?^Xp + 1/32 i?^X,] + 2/3 B(xp, Xpi 0) i?^ Xp 

+ 5/3 B(xp, Xp; 0) LI i?^ xp + B{xp, Xpl 0) [-2 XixZoXp - 4 Xi^?^ Xp + 4 XiR; Xp + 3 Xi<] 

+ B{xp, Xp; 0) LI [-2/3xZiXp - 4/3 R^^ Xp + 4/3 R^ xl + 1/2 i?^ Xp - 1/6 i?^ X,] 

+ B(xp,xp;0)i;: [4xix^^i + 8xii??^Xp-8xii?pXp] + 4B(xp,xp;0)L'7(i?^)' 

+ B(xp,Xp;0)iS [4/3x^^2 + 8/3 i??^Xp- 8/3 i?^Xp] +5(xp,Xp;0)' [- 1/ 18 Xp + V 18 i?;i?^Xp 

+ 1/288 + l/18B(xp,Xp;0)B(Xp,X,;0) [i?^pi?^Xp-i??3^pXp] 

- i/i8S(xp,Xp;0)B(xi,X3;0)i?^^i?^Xp + i/i8B(xp,xp;0)C(xp,xp,Xp;0) (i?^)'xp 

- l/16B(Xp,Xp;0,£)7ri6 [xii?p+i?pXj + 8 B(Xp, X^; 0) i?^p^i?X46p 

+ 8/3B(xp,Xr,;0)Lgi?^j,^i?X46p - 6 5(xp.,Xp.;0)iSxiXp. - 2B{xps,Xps\^)^lxls 

+ i2B(xp.,Xp.;0)igxiXpa + 4B(xp«,Xp.;0)iSxp. + 2B(x„,x^;0)iSxii?:;i3X^ 

+ 2/3B(x^,x,,;0)Lgi?:;i3X^ - 4B(x^,x^;0)igxi^:;i3Xr, + 4B(x^,x,;0)L;i?|n,,^(i?;46i)2 

- B(x^,x,;0)Lg [4/9i?:;i3xl + 8/9i?;;i3X^] + 1/144 B(xi,Xi;0)B(x3,X3;0)i?fi?f 

- 8B{xi,Xs;0)L'7RiRi - 8/3B{xi,X3;0)LlRiRi + 4C{xp,Xp,Xp;0) L^xiRtxp 

+ 4/3C(xp,Xp,Xp;0)iSi?^Xp - 8C'(xp,Xp,Xp;0)LSxi^pXp - 8/3C'(xp,Xp,Xp;0)i^i?^Xp 

+ H^{1,Xp, Xp, X13; X13) [1/8 Xp " 1/16 Xi3 - 1/36 [RPJ^ xi3 + 1/36 R^^R; xi3 - 1/144 (R^f xis] 

+ H^{1, Xp, Xis, X3a; X13) [-1/12 Rlr, Xqs + 1/24 R% xi3 - 1/16 R% Xp + 1/48 R^^ x, + 1/24 R^ Xps] 

- 1/477^(1, Xp4,X<?6,X46;Xi3)X46 + ^^^(1, X»j, Xp, Xi3; X13) [1/36 i?P^i?^i3 Xi3 - 1/72 i?pi?;;i3Xi3] 

- 1/144 i?^(l,x^,Xr„Xi3;Xi3)(i?;i3)^Xi3 - l/48i7^(l,x„,Xis,X3s;Xi3) [R/,sXp~ R^sXq 

- R^ps Xv - R^i3 Xs] + L[^{l,xi,Xi3,X3',Xi3) [l/36-R^^-R^Xi3 - l/18i?3^-R?^Xi3 

- 1/72 i?ii?3Xl3] - l/8-ff^(l,Xl3,Xl3,Xl3;Xl3)Xl3 - 3/8if^(l,Xl4,X34,X4;Xl3)X4 

+ /f^(2,Xp,Xp,Xi3;Xi3) [l/36i?^^i?^Xi3-l/72i?^i?^Xi3] - 1/72 if^(2, Xp, Xr,, Xis; X13) i?Xi3 Xi3 

+ i?^(2,Xp,Xi3,X9;Xi3) [l/36i?^^i?^Xi3-l/72ii^i?^Xi3] + 1/24^/^(2, Xp, Xi., X3.; X13) iipXp. 

- 1/144 if^(5,Xp,Xp,Xi3;Xi3)(i?^)'xi3 - l/72if^(5,Xi,X3,Xi3;Xi3)i?fi?^Xi3 

+ H'" (1, Xp, Xl., X3a; X13) [-5/24 RP^ XqXis + 1/12 R^r, XpsXia + 1/16 R% XgXis + 7/48 R^^ xisXs 

- 1/24 i?^ XpsXls] + l/4i?^'(l,Xp4,X96,X46;Xl3)Xl3X46 + -H'^' (1, Xrj, Xp, Xl3; Xis) [l/9 i?^^i??3 Xl3 

- 1/36 i?^,i?;i3X?3 + 1/18 i?^piipX?3 + 1/72 iJ^iJ^isxL] + ^f^'(l,X^,Xr„Xi3;Xi3) [1/9TO'X?3 

+ 1/9 i??3i?;;i3X?3 + 5/144 (i?;;i3)2x?3] + i?^'(l,Xr„Xi.,X3.;Xi3) [-l/16i?^,XpXi3 -5/48i?^,XqXi3 

+ l/6ii?,ii;,3X?3-l/48i?:;j,,Xr,Xi3-l/48i?;i3Xi3xJ + ^f^'(l,Xi,Xi3,X3;Xi3) [-l/36i?^,i?^x?3 

+ 1/6 i?^,i??^x?3 + 1/72 i?Ji?^X?3] + ^?'"(l,Xi3,Xp,Xp;Xi3) [-l/8xpXi3-3/16x?3-l/36(i??/x?3 
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+ l/l2Rl^Rlx%-ll^%{RlfxU + l/8i?^'(l,Xi3,Xi3,Xi3;Xi3)x?3 

+ 3/8i/^'(l,xi4,X34,X4;Xi3)xi3X4 + i^"" (2, Xp, Xp, Xi3; Xis) [-l/36i?P^i?^x?3 + 5/72i?^i?^X?3] 

+ i/^'(2,Xp,X„Xi3;Xi3) [1/18 i?2pi?^x?3 + 1/72 i?Xi3X?3] 

+ i/^'(2,Xp,Xi3,X9;Xi3) [-l/36i?«Xx?3 + l/72i?^i?^X?3] " 1/24 ff^'(2, Xp, Xi., X3.; X13) i?p XpsXi3 

+ 5/144 i7^'(5,Xp,Xp,Xi3;Xi3)(i?p)'x?3 + 1/72 i/^'(5, xi, X3, X13; X13) i?fi?3 X?3 

+ i/f '(1, Xp, Xp, X13; X13) [X?3 + 2/9 (Rl^f x?3 - 2/9 Rlr^Rl X?3 + 1/9 X?3] 

+ i/3i?f (i,Xp,xis,X3s;xi3)^^^-RsgpX?3 - i/3^^f (i,Xps,x?3,Xp;xi3)-R?^x?3 

- l/3i?f (l,Xps,X93,X»j;Xl3)-R?3-Rsp^Xl3 + l/9i?f (l,Xl3,Xp,Xr,;Xl3) [^g^^^l3X?3 

- RIj,RIj,^RIxU - <'(i,xi3,x„x,;xi3) [1/9 i?^3i?;;i3x?3 + 1/18 Ki3)'x?3] 

+ i?r(l,Xi3,Xi,X3;Xi3) [l/9i?P,i?^x?3-2/9i?i,i??^X?3] + 1/9 (3 X?3 

~ ^p-Rp X13] + i/9^^f (3,xi3,XpiX9;xi3)-Rp,,^pXi3 + i/9-H^f (3,xi3,Xp;X»7;xi3)-Ri3-Rpq,,-RpXi3 

- 1/18 ijf (7, xi3,Xp,Xp;Xi3) (-Rp)^Xi3 - 3/8i7|| (l,Xp,Xp,Xi3;Xi3)Xi3 

- l/8i/fi'(l,Xp,Xls,X3s;Xl3)i?^^i?:gpX?3 + l/8i?2^i'(l,Xps,Xg.,Xp;Xl3) [^?^X?3 + ^s^X?3-^pX?3] 
+ l/8i?|i (l,Xps,X9s,X<j;Xl3)-Rp,,-RspgXl3 + 1/8^^21 (l,Xps7X?s,X»?;Xl3)^13-Rpq,,-fi'sp>,Xl3 

- l/8i?2''i'(l,X^,Xis,X3s;Xi3)i??3^.%^^3,X?3 + ^^s""!' (1 , Xi3, Xp, Xp; Xis) [3/16 x?3 + 1/12 (i??,)'x?3 

- 1/12 i?^^i?^X?3 + 1/48 (iip'xL] + ^^2''i'(l,Xi3,Xp,X,;Xi3) [-l/12i?P^i?;i3X?3 + l/24i?^i?^i3X?3] 

+ l/48i7i;'(l,Xl3,Xr„X„;Xl3)(i?;;i3)'x?3 + ^^2''l'(l,Xl3,Xl,X3;Xl3) [-l/12i?P^i?^X?3 + l/6i?3,,i??^X?3 

+ 1/24 x?3] + 3/8i7i'i'(l,Xi3,Xi3,Xi3;Xi3)x?3 + 9/8 Hfi'(l, X4, Xi4, X34; X13) X?3 

+ 3/4i7|; (l,X46,Xp4,X<?6;Xl3)Xl3 - 1/8 (3, Xps,Xp,X9s;Xl3)^p Xl3 

+ i/2''i'(3,Xi3,Xp,Xp;Xi3) [-l/12i?^„i?^x?3 + l/24i?=i?^X?3] + ^^2''i' (3, Xi3, Xp, X,; X13) [-1/12 x?3 

+ l/24i?^i?^X?3] + l/24i?2''l'(3,Xl3,Xp,Xr,;Xl3)i?Xl3X?3 + l/48i?2''l'(7,Xl3,Xp,Xp;Xl3)(i?p)'x?3 

+ l/24i7f;(7,xi3,Xi,X3;Xi3)i?fi?^X?3- (33) 



IV. DISCUSSION AND CONCLUSIONS 

The analytical formulas for the NNLO contributions 
to the decay constants are very complicated, especially 
when the valence quark masses are nondegenerate. The 
practical usefulness of those results thus depends com- 
pletely on the availability of convenient numerical imple- 
mentations. In the near future, we will make such an 
implementation available [l4j . 

The results, as given in the preceding section, depend 
on nine 0{p^) and five LEC:s. Eventually, their 

proper values should be determined from lattice QCD by 
a fit of the NLO -I- NNLO formulas to the simulation 
data. At the present time, these input parameters have 
been taken from the continuum work in NNLO xPT of 
Ref. [31 . The combination of LEC:s used in this paper 
corresponds to 'fit 10' of Ref. 15J, which had Fq — 87.7 
MeV and ^ = 770 MeV. The parameters which were not 
determined by that fit, namely Lp, L4 and Lg, have been 
set to zero, as have all the K^. This has been done at a 
scale of /J, = 770 MeV. It should be noted that Lq cannot 
be determined from experiment, but is obtainable from 



partially quenched simulations. Some recent results on 
L4 and Lg may be found in Ref. [T^ . 

The graphical presentation of multidimensional func- 
tions is a more immediate problem. In the most general 
case, the decay constant of a charged pseudoscalar me- 
son in PQxPT is a function of two valence and three sea 
quark masses. The plots presented in this section should 
therefore be viewed as an attempt to cover the potentially 
most interesting parts of the parameter space. They also 
serve as a consistency check for the formulas in the pre- 
ceding section, as the cases with fewer different masses 
should be obtained numerically as limits of the more gen- 
eral cases. 



A. Numerical Results for dsca ~ 1 

To this end, the result for dvai — 1 and dsca = 1 has 
been plotted along different rays in the Xi"X4 plane, char- 
acterized by an angle 9 which is defined according to 

tan(0) = X4/X1, (34) 
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such that 9 = 45° corresponds to the unquenched case 
of equal sea and valence quark masses. The typical sit- 
uation in lattice QCD simulations, where sea quarks are 
heavier than their valence counterparts, is attained for 
9 > 45°. The shift in the decay constant due to the com- 
bined NLO + NNLO contributions has been plotted in 
Fig. 121 for different values of 9. The quantity plotted, A, 
denotes the relative change in the decay constant and is 
defined in accordance with Eq. (|24ll , as 



phys 



= Fo (1 + A) , 



(35) 



Thus, when plotted along a curve of the form (|34|1 . A 
is expected to vanish in the chiral limit, which is indeed 
borne out in all plots presented in this section. This 
shows that the NNLO shifts to the decay constant in 
PQxPT do not produce unphysical logarithms. On the 
other hand, the curve labeled 'A' in Fig. [21 is plotted for 
a constant sea mass with Xi = 0.125 GeV^, and in that 
case A approaches a constant as xi ^ 0. 
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FIG. 2: The combined NLO and NNLO shifts of the decay 
constant, plotted for dvai = 1 and dsoa = 1, for values of Q 
between 15° and 75°. The curve 'A' corresponds to xa = 
0.125 GeV^. In order to reduce clutter, most of the lines in 
this plot, with the exception of the 15° line, have not been 
drawn all the way to the origin. 
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FIG. 3: The NLO shift of the decay constant, plotted for 
dvai = 1 and dsca = 1, for values of Q between 15° and 75°. 
The curve 'A' corresponds to X4 = 0.125 GeV^. 



+ 

C\J o 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
Xi [GeV^] 

FIG. 4: The combined NLO and NNLO shifts of the decay 
constant, plotted for dvai = 1, dsea = 2, 6 = 60°, and a strange 
sea quark mass parameter z between 1/2 and 10. 



The results along the lines H34(l at NLO is shown in 
Fig. 121 The chiral expansion for the decay constants con- 
verges better than for the masses as can be seen by com- 
parison of Figs. 121 and with the corresponding ones in 
Ref. 6]. Note, however, that the dependence of A on the 
numerical values of the LEC:s is strong and can change 
the behavior of the NNLO results quite considerably. It 
should also be kept in mind that many of the curves 
plotted in Fig. (21 reach far beyond the expected radius of 
convergence of PQxPT. 



B. Numerical Results for dsca = 2 

The physically more interesting case of dgoa — 2, where 
the strange sea quark mass may deviate from that of the 
u, d quarks, can then be accounted for by the introduc- 
tion of an additional parameter 2:, which is defined as 

^ = X6/X4- (36) 
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The change in A when the strange sea quark mass is 
allowed to vary is shown in Fig. 21 In general, the NNLO 
effects appear to become larger as the strange sea quark 
mass is increased. Fig.0]also shows that Eq. (|28|) reduces 
numerically to the c?vai = 1, ^sca — 1 result of Eq. H27(l 
when all sea quark masses become equal, i.e. z —^ I. 



0.5 




Zi [Gev2] 



FIG. 5: The combined NLO and NNLO shifts of the decay 
constant, plotted for dvai ~ 2, dsoa = 2, z — 5 and 9 = 60°. 
The strange valence quark mass parameter x varies between 
1 and 2. 

It is also instructive to investigate the NNLO correc- 
tion to the decay constant for dvai = 2 and dgca = 2. 
The behavior of A when the strange valence quark mass 
deviates from that of the u, d quarks is shown in Fig. \^ 
where the results have been parameterized in terms of 

X = Xs/Xi- (37) 

It is evident that A changes rapidly with increasing x for 
values of xi above about 0.15 GeV^. Also, Eq. is 
seen to reduce correctly to the dvai = 1, rfsca = 2 result 
when X ^ 1. 



C. Numerical Results for d^ca ~ 3 

In the case of dgoa = 3, and Xn are non-trivially 
related to the sea quark masses, and consequently the 
sea quark sector is more complicated. The variation in 
the mass of the d quark with respect to the u quark, in 
the sea quark sector, has been parameterized by 

y = X5/X4- (38) 

Sample plots of A for dvai — 1, with all three sea quark 
masses different, are shown in Fig. 1^1 Those plots also 
demonstrate numerically the consistency of Eq. (|29|l with 
the dsca = 1,2 results. 




-0.02 - 



0.05 0.1 0.15 0.2 

Xi [Gev2] 

FIG. 6: The combined NLO and NNLO shifts of the decay 
constant, plotted for dvai ~ 1, dsca = 3, 6 — 60° and z — 5. 
The mass parameter y of the down quark in the sea sector 
ranges between 1 and 2. 



D. Conclusions 

A generic feature of the plots presented in this section 
is that the curves for A show a pronounced dip around 
Xi =0.1 GeV^. This indicates a strong cancellation be- 
tween the NLO and NNLO contributions to A. It should 
be noted that such a feature is apparently not exhibited 
by the analogous expressions for the pseudoscalar meson 
mass 0. This cancellation also depends strongly on the 
choice of the LEC:s. Many more plots can of course be 
produced but those presented give a first indication of the 
size of the corrections and how they vary with the differ- 
ent quark masses used as input. We have not attempted 
any fit of our results to the available lattice data; Many 
simulations are performed with only two flavors of sea 
quarks, and in addition we need extrapolations to zero 
lattice spacing and infinite volume at each quark mass to 
apply the present formulas. Work for the case with two 
sea quark fiavors is in progress. 

In conclusion, we have calculated the decay constants 
of the charged, or off-diagonal, pseudoscalar mesons to 
NNLO in PQxPT and presented analytical as well as nu- 
merical results for a variety of different combinations of 
quark masses. The NNLO contributions were found, as 
expected from previous work in NNLO xPT [ij, llSl , 
to be sizable even though there is a tendency toward can- 
cellation with the NLO result. As the results depend on 
a number of largely unknown LEC:s, statements about 
the convergence of the chiral expansion have to be post- 
poned at this time. The NNLO effects are definitely non- 
negligible at presently used quark masses in Lattice QCD 
simulations. 
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